Particle physics catalysis of thermal Big Bang Nucleosynthesis 
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We point out that the existence of metastable, r > 10^ s, negatively charged electroweak-scale 
particles {X~) alters the predictions for lithium and other primordial elemental abundances for 
A > 4 via the formation of bound states with nuclei during BBN. In particular, we show that the 
bound states of X~ with helium, formed at temperatures of about T — 10*K, lead to the catalytic 
enhancement of ''Li production, which is eight orders of magnitude more efficient than the standard 
channel. In particle physics models where subsequent decay of X~ does not lead to large non- 
thermal BBN effects, this directly translates to the level of sensitivity to the number density of 
long-lived X~ particles (r > 10^ s) relative to entropy of nx~/s ^ 3 x 10~^^, which is one of the 
most stringent probes of electroweak scale remnants known to date. 



Standard Big Bang Nucleosynthesis (SBBN) is a well- 
established theory that makes predictions for elemental 
abundances of light elements, H, D, He and Li, as func- 
tions of only one free parameter, the ratio of baryon to 
photon number densities. Agreement of the observed 
abundances for D and ^He with the SBBN predictions 
that use an additional CMB-derived [l| input value of 
rib/s = 0.9 X 10^^° serves as a sensitive probe of New 
Physics. Indeed, over the years a finessed SBBN ap- 
proach has led to constraints on several non-standard sce- 
narios, including constraints on the total number of ther- 
mally excited relativistic degrees of freedom P| . This first 
class of constraints results essentially from the destor- 
tion of the time- frame for the n/p freeze-out, affecting 
mostly the abundance of ^He. The second class of con- 
straints comes from the nonthermal BBN, that results 
from the late decay of metastable heavy particles 0] . If 
such decays happen during or after BBN, they trigger 
electromagnetic cascades that affect light elements via 
photodissociation, or lead to non-thermal nuclear reac- 
tions by fast hadrons produced in the decay, resulting in 
significant modifications of D, ^Li ''Li and '^He/D abun- 
dances. The BBN constraints have been instrumental in 
limiting some variants of supersymmetric (SUSY) mod- 
els with long-lived unstable particles. For example, late 
decays put constraints on the energy density of unstable 
gravitinos, limitirig it to be less than ~ 10~^^ relative to 
entropy xGeV 0,0,01 for some selected range of lifetimes 
and SM branching ratios. It is important to establish 
whether the late decay of heavy particles can possibly 
"cure" [3, 0, @| what is known as the lithium problem, a 
statistically significant and persistent discrepancy of the 
SBBN prediction for ''Li from the roughly twice smaller 
observational value over a wide range of metallicities. 

The purpose of this Letter is to show that in ad- 
dition to the change of the timing for BBN reaction 
and non-thermal processes, there is a third way particle 
physics can significantly affect the prediction for primor- 
dial abundances of light elements. This new way consists 
of the catalysis of thermal nuclear reactions by heavy 



relic particles that have long-range (electromagnetic or 
strong) interactions with nuclei. In particular, we show 
that the cosmological presence of metastable charged par- 
ticles, called X~ hereafter enables the catalyzed BBN 
(CBBN) via the formation of bound states between light 
nuclei and negatively charged particles X~ . These bound 
states form in the range of temperatures from 1 to 30 
KeV, changing the standard nuclear reaction rates but 
more importantly opening new channels for thermal re- 
actions and changing the abundance of Li and other ele- 
ments with A > 4. The most significant difference is seen 
in the ^Li production mechanism. 



SBBN : *He D ^ ^Li + 7; 
CBBN: (4HcX")+D 



^Li 



Q = 1.47MeV (1) 
X-; Q ~ LlSMeV, (2) 



and, as we are going to show, the cross section for the 
CBBN channel is enhanced by eight orders of magnitude 
relative to SBBN. In ^ and below, {NX) denotes the 
bound state of a nucleus N and X~ . In the remainder 
of this Letter, after a brief review of the relevant bound 
states (NX), we analyze the cross section and thermal 
rate for reaction ([2]), make a prediction for ^Li in CBBN, 
put constraints on some of the particle physics scenarios, 
and point out new CBBN mechanisms for ''Li depletion. 

Properties of the bound states. Here we assume that 
the electromagnetic force is acting between nuclei and 
X~ particles, and calculate properties of the ground 
states using the variational approach. Binding energies 
Eb, average distances between the center of the nucleus 
and X~, and the "photo-dissociation decoupling" tem- 
peratures To, are summarized in Table I where the con- 
straint mx- S> TTiM is also imposed. The main uncer- 
tainty in Table I comes from the charge distribution in- 
side the nucleus as the naive Bohr orbit oq = {Zamf^)~^ 
can be well within the nuclear radius. It leads to a reduc- 
tion of the bound state energies relative to the Bohr-like 
formula, E'^ = from ~ 13% in (^HeX) to 50% 

in (^BeX). Realistic binding energies are calculated for 
two types of nuclear radii assuming a uniform charge dis- 
tribution: for the simplest scaling formula = 1.22^3 , 
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larger than the Hubble rate. For the same reason the 
abundance of bound states with rare light elements is 
very low, nj^x- /n-x- ^ 10^^, where N ^■^He, D and T. 
As we are going to show, an actual abundance of X~ at 
T ~ 8 KeV has to be less than 10 which makes the 
overall impact of these bounds states on BBN negligible. 




TABLE I: Properties of the bound states: Bohr ao and nuclear 
radii Rn in fm; binding energies Eb and "photo-dissociation 
decoupling" temperatures To in KeV. 

and for the nuclear radius determined via the the root 
mean square charge radius, Rnc = (5/3) ^/"^-Rc with ex- 
perimental input for Rc where available. Finally, as an 
indication of the temperature at which (NX) are no 
longer ionized, we include a scale Tq where the photo- 
dissociation rate rph(r) becomes smaller than the Hub- 
ble rate, rph(To) = H{Tq). It is remarkable that sta- 
ble bound states of (^BeX) exist, opening up a path to 
synthesize heavier elements such as carbon, which is not 
produced in SBBN. In addition to atomic states, there 
exist molecular bound states (NXX). The binding en- 
ergy of such molecules relative to (NX) are not small 
(e.g. about 300 KeV for ("HeX-X")). Such neutral 
molecules, along with (^BeX) and (^BeXX), are an im- 
portant path for the synthesis of heavier elements in 
CBBN. One can easily generalize Table 1 for the case of 
doubly-charged particles, which was recently discussed in 
Q in connection with the dark matter problem. 
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FIG. 1: Fraction of X~ locked in the bound state with *He: A. 
Realistic result based on Boltzmann equation, B. Saha-type pre- 
diction with a rapid switch from to 1 at T ~ 8.3 KeV 

The initial abundance of X^ particles relative to 
baryons, Yx(t <^ t) = nx-fni,, along with their lifetime 
T are the input parameters of CBBN, and it is safe to 
assume that Yx <^ 1. The most important catalytic en- 
hancement results from the bound states ("'HeX). Their 
abundance relative to the total abundance of X^, Ybs = 
n-Bs/fT-x- is calculated using the Boltzmann equation 

dY 

- HT—^ = {(7recv)nB_c " (o-p/i?;) VbS"-7 (3) 

along with the photoionization and recombination cross 
section. The result is shown in Fig. 1, where a significant 
deviation from the naive Saha equation is observed as the 
recombination rate of X~ and ^He is only marginally 



FIG. 2: SBBN and CBBN mechanisms for producing ''Li. 

Photonless production of ^Li. The standard mecha- 
nism for ^Li production in SBBN is "accidentally" sup- 
pressed. The D-'*He cluster description gives a good 
approximation to this process, and the reaction rate 
of ([1]) is dominated by the E2 amplitude because the 
El amplitude nearly vanishes due to an (almost) iden- 
tical charge to mass ratio for D and ''He. In the E2 
transition, the quadrupole moment of D-^He interacts 
with the gradient of the external electromagnetic field, 
Vint = Qij'^iEj. Consequently, the cross section at BBN 
energies scales as the inverse fifth power of photon wave- 
length A — Lo^^ ^ 130 fm, which is significantly larger 
than the nuclear distances that saturate the matrix ele- 
ment of Qij , leading to strong suppression of ([T]) relative 
to other BBN cross sections [9]. For the CBBN process 
([2|) the real photon in the final state is replaced by a 
virtual photon with a characteristic wavelength on the 
order of the Bohr radius in (^HeX^). Correspondingly, 
one expects the enhancement factor in the ratio of CBBN 
to SBBN cross sections to scale as (agw)"^ ^ 5x10^. Fig- 
ure 1 presents a schematic depiction of both processes. 
It is helpful that in the limit of i?jv <S ooj we can ap- 
ply factorization, calculate the effective ViEj created by 
X~, and relate SBBN and CBBN cross sections with- 
out explicitly calculating the (D^HejQy |^Li) matrix el- 
ement. A straightforward quantum- mechanical calcula- 
tion with ViEj averaged over the Hydrogen-like initial 
state of (■*HeX) and the plane wave of ^Li in the final 
state leads to the following relation between the astro- 
physical S'-factors at low energy: 

c c ^ Pf"'0 1^1 , ™D V 

Here ao is the Bohr radius of ('^HeX^), pf — 
(2TO6Li(QcBBN + E)y^^ is the momentum of the out- 
going ^Li in the CBBN reaction, and ui is the photon 
energy in the SBBN process, uj — Qsbbn + E. For 
E <^ Q the value of the final momentum of the ^Li nu- 
cleus is Pf ~ (1.8fm)~^. Throughout the whole paper, 
c = fi, = 1. The S factor is defined in the standard way, 
by removing the Gamow factor G from the cross sec- 
tion: S{E) = Ea/G. A somewhat more sophisticated 
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approach with the use of an optimized wavefunction for 
{^KeX) in the initial state and the Coulomb wavefunc- 
tion in the final state, with the radial integrals calculated 
with the constraint jr^ — reLij > 2 fm, renders the result, 
which is very close to a naive estimate: 

^cbbn(O) ^ 6 X 10^ X c 0.3 MeVbn. 

(5) 

Here the SBBN value of 5sbbn ^ 1-8 x IQ-^eVbn is 
taken at some matching scale Eq — 0.41 MeV [1], which 
is sufficiently close to 0, yet is in the regime where SBBN 
cross sections can be reliably measured. The result of the 
factorized approach ([5]) should be within a factor of a few 
from the exact answer, with the largest errors presumably 
associated with the Rn/qq expansion. Using the Gamow 
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FIG. 3: Logio(^Li) plotted as a function of T in KeV for different 
cfioices of T and Yx- A. r = oo and Yx = 10~^, B. r = 4 X 10^ s 
and Yx = 10"^, C. r = oo and Yx = 10'^, D. t = 4 x 10^ s and 
Yx = 10-^ 

energy that accounts for the charge screening in {^HcX ) 
and the change of the reduced mass, 



e: 



Gamow 
SBBN 



5249 KeV 



e: 



Gamow 
CBBN 



1973 KeV, (6) 



we calculate the thermally averaged cross section for the 
CBBN processes at small temperatures T ~ 10 KeV em- 
ploying a saddle point approximation: 

(crew) ~ 2bn X T-2/3exp(-23.7T-i/3) 

= 1.8 X 10^ X rg"^/^exp(-5.37Tg"^/^). (7) 

In the second line, the cross section is converted to cus- 
tomary units of N^^ cm^s~^ , and Tg is temperature 
in units of 10^ K. CBBN rate (O is to be compared with 
the SBBN expression 



asv) ~ SOTg ^/%xp(-7.435rg 



-l/3x 



(8) 



and the enhancement of eight orders of magnitude for 
CBBN is traced back directly to ([S]). Although tremen- 
dously enhanced relative to (O, the CBBN rate ([7]) is by 
no means larger than the usual photonless reaction rates 
known to SBBN. 

Catalytic enhancement o/ ^Li at 6 < T < 12 KeV. 
Armed with the rate ([7|), we write down the evolution 
equation for ^Li abundance at T < 12 KeV: 



In this formula, D and ®Li are the Hydrogen-normalized 
abundances of these elements, Up and H are the 
temperature-dependent concentration of free protons 
and the Hubble rate, and (cTpv) is the rate for 
^Li+p ^■^He-l-^He responsible for the destruction of ^Li 
[i. We use the output of the SBBN code for D, "He as 
function of T and solve the CBBN subset of equations ([3]) 
and ([9]) numerically. In fact, since n^siT) significantly 
differs from zero only below 9 KeV, the variation of '''He 
and D abundances is negligibly small compared to the 
freeze-out values (D= 2.4 x 10~^, Yp = 0.25). Several of 
the numerical solutions to this equation for different in- 
put values of Yx and r are plotted in Fig. 2. The increase 
in time of ^Li below 10 KeV is a direct consequence of 
CBBN. One can clearly see that the with Yx ^ 0(0.01) 
one could convert up to a few % of D to ^Li. Taking the 
limit of large lifetime (r > lO^s), and using the limiting 



abundance of primordial Li of 2 x 10 
remarkable sensitivity to X~: 



\\, we get a 



Yx < 3 x 10" 



nx-/s < 2.5 X 10 



-17 



(10) 



A scan over the lifetime parameter produces the ex- 
clusion boundary on the (r, logj^gYx) plane plotted in 
Fig. 3. Also shown is the natural range for Yx cal- 
culated using the standard freeze-out formula with the 
input annihilation cross section in the range of aaxmV = 
(1. - 100.) X 2Tra'^m^^ and mx = (0.1 - l)TeV. Taken 
literally, the ^Li overproduction constrains the heavy 
charged particles to the range of lifetimes 



r < 5 X lO^s, 



(11) 



unless (TunnV is tuned to larger values. In fact, both re- 
sults, (fTO|) and (fTTjl . should be interpreted as hmits of 
sensitivity to {t,Yx) via the CBBN enhancement. It is 
true that excessive energy injection at later time may re- 
duce the amount of ^Li generated in CBBN, and as a 



consequence relax the constraints of Eq. (jlOp [llj . The 
combined analysis of two effects, the CBBN enhancement 
of ^Li and impact of energy release on elemental abun- 
dances, has to be done by combining the result of this 
work and previous papers on unstable particles in BBN 
[1, 0, [1] . However, such analysis is necessarily model de- 
pendent, and falls outside the scope of the present Letter. 
In many models, however, the decay of X~ is not involv- 
ing huge energy release, and then both pUj) and (fTTI) hold 
unconditionally. 

CBBN points towards an intriguing possibility for al- 
leviating the lithium problem. As is evident from Fig. 3, 
lifetimes of X~ between 2000 and 5000 seconds can cre- 
ate a primordial source for ^Li thus removing some ten- 
sion between observations of ^Li in low metallicity sys- 
tems and explanation of its production using cosmic rays 
But more importantly, such lifetimes of X~ can 
lead to an overall reduction of ^Li abundance. At the 
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WMAP-determined baryon density, most of ''Li is pro- 
D(nBs(CTcw) + nHe(crsu}) Li np{apv){%) duced via ^Be. At T - 30 KeV and r ~ (2 - 5) x 10^ s 
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X~ is still present in sufficient quantities to trigger the 
formation of (^BeX). Once formed, (JJieX) facilitates 
several reaction rates that destroy ^Be. The decay of 
X~ within the bound state leads to further depletion 
of ''Be and ultimately of ^Li. Moreover, if the decay of 
X^occurs via a virtual VF-boson, there exists a fast elec- 
troweak capture process of X~ on "^Be, which also leads 
to an overall depletion of ^Li. 

Usually, a serious obstacle on the way to solve Li prob- 
lem is the abundance of D and ^He/D ratio 0, H, 0| 
affected through large energy release, but we emphasize 
that CBNN mechanism itself is inconsequential for lighter 
elements such as D or ^He as the abundance of any bound 
states of X~relative to baryons is not larger than 10^^ 
once the ^Li bounds are satisfied. 




FIG. 4: Constraint on initial abundance of X~ . logi^iYx) is 
plotted against t in units of 10^ s. The thick horizontal line is the 
asymptotic exclusion boundary JTOJ. 



Prior to this work, models of type II were believed 
to be of no consequence for BBN. Type II models are 
not involving significant energy release, and therefore 
our results based on CBBN enhancement, Eqs. (fTO|) 
and pT|) . become strict limits rather than levels of sen- 
sitivity. Examples of type II include an important case 
of near-degenerate states of lightest (LSP) and next-to- 
lightest (NLSP) SUSY particles, neutralino-chargino or 
stau-neutralino. Using recent calculation Ref. of 
stau lifetime, we can translate (jlip into the limit on 
minimal mass splitting in the stau-neutralino system, 
ruf — > 70 MeV. Along the same lines, one can 
derive constraints on the minimal mass splitting between 
neutral and charged particles at the first excited Kaluza- 
Klein level in models with "universal" extra dimensions. 

To conclude, thermal catalysis of BBN reactions repre- 
sents the novel idea of how particle physics can alter the 
predictions for light elemental abundances. The cataly- 
sis of ^Li via the production of (*HcX^) bound states is 
explicitly demonstrated, and the enormous enhancement 
in the catalyzed channel furnishes the sensitivity to X~ 
normalized on entropy at 3 x 10~^^ level. I would like to 
thank C. Bird and K. Koopmans for assistance in some 
calculations and useful discussions, and acknowledge use- 
ful conversations with Drs. R. AUahverdi, R. Cyburt, K. 
Olive, A. Ritz and Y. Santoso. 



Particle physics applications. There are two generic 
possibilities for a long-lived X^: either the decay of 
X~ is suppressed by tiny couplings, or it decays to an- 
other neutral relic state X^ with small energy release. 
Schematically, these possibilities can be represented by 

Type I : X" SM"[X°], AE Mx, (12) 
Type II : X~ X° + [ly] ; AE < few MeV, 

where the first type corresponds to a decay into charged 
SM state(s) with or without any neutral non-SM relics, 
and the second type is X~ decaying to another neutral 
relic plus an electron and possibly some neutrino states 
with a small energy release. SUSY models have scenarios 
of both types. An example of type I is the charged slepton 
{e.g. stau) decaying into gravitinos or decaying due to 
i?-parity violation without leaving any SUSY remnants. 
Since the amount of energy released in the decay oiX~ of 
type I is large, we cannot restrict this scenario only to the 
CBBN-enhanced output of ^Li as non-thermal BBN may 
provide equally important changes in the elemental abun- 
dance. In case of very late decays the constraints coming 
from e.g. diffuse 7-ray and microwave backgrounds can 
become important. 
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